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Abstract 

This paper describes a modelling technique used to explore three dimensional (3D) image ^radiance distributions formed 
by high numerical aperture (NA > 0.5) lenses in homogeneous, linear films. This work uses a 3D modelling approach that 
is based on a plane-wave decomposition in the exit pupil. Each plane wave component is weighted by factors due to 
polarization, aberration, and input amplitude and phase terms. This is combined with a modified thin-film matrix technique 
to derive die total field amplitude at each point in a film by a coherent vector sum over all plane waves. Then the total 
inadiance is calculated. The model is used to show how asymmetries present in the polarized image change with the in- 
fluence of a thin film through varying degrees of focus. 

1. Introduction 

The demand for highly integrated electronic devices has motivated investigation into better lens resolution in micro- 
photolithography. Lens resolution is roughly proportional to k/NA; therefore, we desire to maximize NA and minimize k. 
Photoresist, source, and lens materials issues limit the operating wavelength much below 248nm. NAs are design and 
budget limited, with a possible practical maximum of0'95. In this paper, we will show that sophisticated models are needed 
to understand fundamentals of imaging in the high NA regime. 

Traditional models used in micro-photolithography are based on scalar image formation under the Fresnel 
approximations 1 . The central approximation is to treat spherical lens and pupil surfaces as pure quadratic paraboloids. This 
holds in the low NA regime but it breaks down when the exit pupil diameter is of the same order as the distance from pupil 
to image, i.e., for NA > 0.5. Furthermore, the traditional scalar models cannot ueat polarization. They, in effect, assume 
that each plane-wave component has the same polarization amplitude, with polarization vectors parallel to one another and 
perpendicular to the axis. In a high NA system, there is significant variation of these vectors across the exit pupil. Each 
x, y, and z co m ponent has a different polarization contribution. Furthermore, these scalar models treat irradiance 

within the photoresist film as if the aerial image could be propagated through the film stack as a normally incident plane 
wave. This approximation is reasonable for NAs below 0.5 or 0.6, but at higher NAs, the propagation angles of the electric 
field (E) become significant 

One of the early experiments of Wiener 2 shows that, at oblique angles, light polarized parallel to the plane of incidence ( P 
or TM polarization) exhibits reduced contrast in photographic film compared to light polarized perpendicular to the pt-ine 
of incidence (S or TE polarization). For example, a pair of P-poiarized plane waves propagating within a film at r45° 
produce no interference fringes because their E vectors are perpendicular to each other. Moreover, standing wave elfccts 
in the resist film stack are also different in the two polarizations. This is illustrated in Figure 1 with a plane wave incident 
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on a photoresist film over silicon at 33°, 44°, 58°, and 72° (or 0.55, 0.7, 0.85, and 0.95 NA respectively). The total power 
absorbed in the film, which is proportional to photoresist exposure, is plotted as a function of coated photoresist thickness 
in Figure 1. Although the differences between the polarizations are slight at 0.55 NA, they become appreciable high NAs. 
It becomes apparent (hat any imaging models for high NA must take into consideration the vector nature of light. 


This work presents a 3D vector imaging model for homogeneous, linear, thin films. Previous work 5 concentr ated on a 2D 
approach. Other authors 4 ' 1 have looked at computationally intensive techniques for inhomogeneous films or required 
axially symmetric systems to obtain analytic solutions. In comparison, by limiting this work to homogeneous films, a 
technique can be presented that is easily visualized and computationally rapid. Since the final computation is numeric, no 
restriction is placed on the system symmetry. The model uses the Debye approach 9 to characterize the image field as a 
plane wave-decomposition. This restricts the propagation direction of the plane waves to the cone of rays subtended by the 
exit pupil of a high NA lens with its vertex at the geometrical focus. If a film stack is located at or near the focus, the 
amplitude and phase of each plane wave can be used as input into a thin-film matrix routine to calculate the local field in 
the top film. The total field is then the summation over all plane waves. 

2. Model 

Consider the Kohler projection system shown in Figure 2. It is assumed that the system is lossless, and has a Strehl ratio 
above 0.8. The film of interest is the 1st film, which is homogeneous and linear with thickness dj. All optical elements 
from the source to the entrance pupil have low NA and can be modelled by traditional scalar methods. Additionally, the 
polarization remains constant throughout the optical system and is only altered at the exit pupil. These assumptions allow 
for a direct mapping of the amplitude and phase of the entrance pupil to the exit pupil of the imaging lens. In Figure 2, 
the source is initially polarized along the y axis. If a ray is considered to emerge from the exit pupil, it will now have 
polarization amplitudes for each cartesian component that depend on the exit pupil location. 

The main assumption of this model is that each point in the exit pupil gives rise to a Huygens secondary spherical wavelet. 
If the image fieldTsze is much less than the pupQ diameter and is located abtxinhe axis, W wi^ttscah be represented 
as plane waves propagating toward the the geometrical focus. Since each plane wave from an exit pupil point is normal 
to the localized surface about that point, the plane wave directions are limited to the cone (for circularly symmetric systems) 
formed by the vertex at the focus and the radius of the pupil. This is the Debye diffraction approach*. In the cases con- 
sidered here, a local field size less than 10pm about the axis is sufficient to encompass all object points contributing sig- 
nificant amplitude to features of interest 
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In Figure 3 we illustrate a general propagation vector, k, for a plane-wave as 

“4 AAA* 

k - ^q( cos0sin0x + sin0sin0y+cos0z) 

-fyl(ax + 0y + Yz) , 

where io is the free space propagation constant , and q is the defined as complex index of refraction given by 

q = n — iK , 

where n is the real part of die refractive index and k is the extinction coefficient 


■4 -♦ 

Since 0 is defined as the angle thru k makes with die z axis and 0 as the angle that the projection of k on the x-y plane 
makes with the x axis, the direction cosines have the standard definitions: 

a = cos 0 sin 0 

P = sin 0 sin 0 


y= cos 0 
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Y = >/l-a 2 -p 2 . 


Propa ga tion of a plane wave to the 1st film is illustrated by Figure 4. The pupil is given as a surface in direction cosine 
sp are over the variables (a, (}). Since z - 0 at the center of curvature of the exit pupil, the direction of each piane-wave 
emulating from the exit pupil location (<*>, (Jo) is completely defined by (oto. (Jo)- A coordinate z' is defined as 

z' = z - zq . [2] 

where zo is the distance from the geometrical focus of the system to the top surface of the film stack referenced to the in- 
cident medium. This places the top film surface at z = 0 and the geometrical focus in the incident medium at z = - : 0 or 
7 = 0. If no aberrations are present in the optical system, all plane waves from the pupil are in phase at z = 0. Each plane 
wave will have a phase of — AodoYoZo at z — zo, where subscripts on txo, (Jo, Yo» and T|o reference th$ incident medium. Also, 
the position vector, f , has its origin atjhe top film surface, and the propagation vector in air, ko, is coplan ar with the re- 
fracted propagation vector in the film, 

Each monochromatic plane wave propagating from the pupil to the top surface of the film has the form (with the periodic 
lime factor e-“ suppressed) 

£<a.(S) = £,<a.|S)«-* •' = £*a.|»V% '-*»** ?<«.» . U1 

and 

kQ.? = Mo( a o x+ Po? + Tfo*') • 

a, (5) describes the amplitude and phase on the entrance pupil due to the object diffraction. It is a scalar term due to the 
initial assumption of low NA oo the object side. The ^ term is arises from the requirement that the system has negligible 
energy losses and is aplanatic. It is derived 1 using the conservation of energy law, and it results in a complex amplitude 
at the exit pupil of &( a, |J) -^b . The next exponential term describes any residual aberration in terms of a wavefront phase 
error, W(a, jj). Following this is the phase term due to the location of the top surface of the film. P( a, (1) is the polarization 
vector amplitude distribution across the exit pupiL It can be derived in terms of the propagation vector 10 , ko, since at all 
times it is assumed that the polarization vector is perpendicular to the propagation vector. Prior to the imaging lens, the 
i nitial polarization vector of the electric field is in some state of polarization that can be treated in terms of a orthogonal 
decomposition in x and y. Therefore, one only needs to examine the incoming component polarization vectors along the 
x or y axis. After passing through the lens, a given ray will have its polarization vector rotated accordingly. The cartesian 
components are found by applying the appropriate coordinate transformations. Upon interaction with the thin film surface, 
a plane of incidence about the z axis can be defined by the angle *. This allows the components of the polarization sector 
to be further reduced in terms of contributions due to a component perpendicular to the plane of incidence, S or TE, and 
parallel to this plane, P or TM Table I shows the values of the polarization amplitudes for each component. Note that 
PI + pz + pi - i and that the z component is only comprised of P polarization. 

Figure 5 illustrates a plane wave from the pupil with amplitude and phase E<Jsx, ()), described by Equation 3, arriving at the 
top surface of the film, interface L The total field at any x,y,z point within the top film is found by summing the downward 
and upward fields at dim The derivation presented here departs from the typical thin film-methods by presenting 

this in terms of the fields at bottom of the film, interface II, and then relating the result to the incident field at interface I 
through matrix formalism. The incident and reflected fields at interface II give for the field in the 1st film. 




Downward Field, 
Upward Field, 
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where z" = d x - z. 


The following relationships are found from Snell's law in direction cosine notation: 


a i =T lo a o 
TllPl=Tlo3o ’ 


and 


»ll Yi = 



[ 7 ] 


The total field is given as 

£j(ot, P) = E* (a, P) + £' (a, p) 

= ( Ej,(a, P) + E HI (a, P) c-** n ) , 1 1 

where 

<KO = koz“^r\]-r\l sin 2 ( cos" ‘y 0 ) • [9] 

The amplitudes in Equation 8 are related to the incident field by using thin-film matrix techniques to derive reflection and 
transmission coefficients for the full film stack and a partial film stack. Since the basics of this technique 1141 are well 
known, only the novel mathematics pertaining to the current discussion will be presented. It is assumed that the index of 
refraction above the first surface, ifo. is real, and that the film in question is the first film layer on an arbitrary film stack. 
Using the notation in Hecht u , a matrix can be defined for each film, j, by 


M; 


COS&; 


i sinSy 
iYj sin cos 


[ 10 ] 


where Y t is defined as the effective index for S and P polarization and 8, is the phase for a film thickness 4. given by 

Yjs^jYj 



5 , = M, 4 Y/ . 


[H] 


A total M matrix can be defined that is the product of all the individual film matrices, as well as a partial M' matrix that 
represents the film stack without the contribution of the 1st film: 


Total Matrix, 


M 


=rb=[ 


*u 

*21 


*12l 
*22 J ’ 


[12] 


and 
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Partial Matrix, M' = rTM =[ m , u "V 2 1 . 

*■ 4 ' I 21 m 22 J 

;-2 


[131 


With respect to Figure 5, the reflection and transmission vector coefficients for the full film stack, r and x respectively, and 
at the 2nd interface, r a and x ff , can be defined with components given by 


and 



^zry 

Ejy 

z ~ E Jt ’ 

[C] 

E 

_ _ trx 

*11 ~ *flx “ *lly ~ u 

c 'illz 

E-tiy 

Ejly 

x - ^ Ul 
X E Jh ’ 

[15] 

E r ]x 

r - r * = r y E Ux 

Erfy 

Euy 

r = ^ 

1 E* ’ 

[16] 

Erllx 

11 = l,x = u > ~ Etfix 

Erlly 

Ejly 

Erl] i 

r "'* ' 

[171 


where the original assumption of film homogeneity and linearity has been invoked to allow the equality of the x and y 
coefficients. 

The coefficients can be expanded following a similar treatment outlined in Hecht in terms of the M and M' matrix elements, 

2T, 


x = 


2]0 _ _ ±U 

F(ynj i + 1 + m 21 + 22 


Klfft'n + Tiiyw'ia - n^2i ~ Tyn a 
^m'tt + Til 

where the dependence of the coefficients on o and p is given by Equation 11. 


T<y»i i + Tq T j/n^ - ff»2i - T/ng 

r T<yni i + ToT/n 12 + «2i + Tyn^ 11 ^m'u + Y\Yjn'n + "*'21 + 22 


18 ] 
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Propagating plane waves in the top film, which are represented by Equations 4 and 5, must also be solutions to the homo- 
geneous wave equation and, therefore. Maxwell's equations. In particular, for a charge free media, it is required that the 
divergence of the field is zero in the incident medium and the film in question, hence 


V . £ * k • E = 0 


Equations 6 and 14 for the incident medium give 

-* n _ T|o <®0 £to + Po n 

*0 • 0 (<*0 E itx + Polity + Yo E Ut ) = x + X, U ’ 


CO] 


Cl] 


which gives 


[ 22 ] 


Ho«>o£» + Po^W) - . 

Equation 20 for the downward field in the film and equation 22 give 

k x ■ E* * t|j (ai E Jlx + Pj E u , y + y, Ejfj ) 

Tlo (®o + Poetry) _ Til Yl 

z fl X !h 

-*TioYo E w niYi^w _ n 

t,t // T //r 

Therefore, 

J*_ T n 0 Yo 

T //l T // TllYl 


[ 23 ] 
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Similarly, for the upward field, 

*1 • £ = nt («l E rllx + Pi £ r//y “ Yl E r!!z ) = r // ^ 1 (°1 E i llx + Pi %y) ” r //z ’ll Yl Etf/z = 0 
By Equation 23, 

ni ( a i Ej/x + Pi Ejly ) * “ Hi Yl £|//z ■ 

Then, by substituting Equation 26 into 25, 

r // n i Yi * - r //z ni Yi fir/* . 

and, 

r Ih ~ ~ r U - 


[25] 

[261 

[27] 

[ 28 ] 


Finally, using Equations 3, 14-17. 24, and 28, Equation 8 can be given for each cartesian component, linking the thin film 
field with each plane wave emanating from the pupil, that is, 

B^a. P) = a, p) Vf e^ (a>p) /yo, p) (e^° + r,,*"* 0 ) , [29] 

£ ry ( a , p ) = **61 <“*■**> e "* 0 ™* 8 £^ ou 3 )>/y /»,« x , P ) ( e ** 0 + r // e ’ < * , " ) ) , [30] 

and 


£ T2 (o, 3) - e* 8 ” e"*’ 17 * 8 ^ p) ^ «*«**•» ^ p) ^ > , [3 1] 

where a, 3, y, and q are referenced to the incident medium. 

Since reflection and transmission coefficients are only defined in terms of either the S or P orientation. Equations 29-31 
are further reduced into S and P components. The total field at any point within die 1st film is then given by the component 
sum of Ej{a, 3) due to the S contributions and the P contributions over the solid angle subtended by the pupil. The nor- 
malized field is given by 
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£ (xou) = K f £ rs <a, P) dfl + K f £ rP (a, P) dQ , 
■'a J n 


[32] 


where K is a the normalization constant 


Since the differential solid angle can be written as 


dQ = 


da dp 


[33] 


the final expanded form for the local electric field components in a film is: 


= K 


II 


£q(a, P) P^q^P) P) ^(a x y -y 1 j _£s_ + r ns e ~ ) da dp 


V7 


II 


£0(0, P) P) e** 0 ^* P) ^(ax+py-y*,,) _V_ (e *K*") + r „ P e~‘* in )dac® 


V7 


x ///» 


[34] 


Ejtxjj) = K 


II 

*a*p 


£p(a, P) fysfa, P) e ^ 0 (a,P) ^^ax+py-rrj) Js_ (e**n + r ns e~ i ** r) ) da dp 
_ /IT *//$ 


1 

•'a *P 


V7 

£o((X, P) fypfot P) g*° ^^iKax+Py-rio) + r[f/>e -‘*X‘" ) )dad$ 

V7 


x ///> 


[35] 


£ r (x^) = K 


II 


£o(a, P) Ptpja, P) e** 0 ^ 8, ^ y -r *) WY . g - «♦(*"> _ j da dp . 


V7 


tW^ltYl 


[36] 


The limits of integration for circularly symmetric systems extend from -NA to NA with the requirement that 


4 


a 2 + p 2 -SM4 . 


[37] 


The localized irradiance at any point within the film can be expressed as, 

I(x,yj) = [£(x,yx) • £*(x,y,z)] /?g<'n 1 ) = ( I Ej^xyj) \ + I Ey(x,yj ) | + |£ t (xo'd)l )£eOli) ■ l J 8] 


This equation is easily evaluated numerically. The minimum sampling size can be determined either by trial and error, 
looking for convergence of the image, or by recognizing that Equations 34-36 represent a Fourier transform. If the image 
function is band limited within the image field, then by sampling theory 14 , the image field sizes, X and Y, are inversely 
proportional to the sampling in the pupil- The minimum direction cosine interval can then be written as 
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[39] 


Aa = 

AP = 


_X 

X 

x_ 

Y 


In the cases of interest in this work the field size is limited to 5\im at \ - 248/i/n giving Aa = 0.05 For a 0.95 NA circ ular ly 
symmetric system this results in 1 100 samples in the pupil, i.e., only 1 100 plane waves are needed to create the image within 
a film. This number is much less than other imaging models based on plane-to-plane propagation. 


3. Results 

Equation 38 was evaluated using a wavelength of X = 248/un and an incident medium of air with the wavefront aberration 
set to W(a, P) = 0. Initially, the image of a point, or the point spread function (PSF), was computed for a non-absorbing 
photoresist film of rfr = 1.8 on a matched substrate for an NA of 0.95. Figure 6 shows the normalized iso-irradiance con- 
tours for meridional slices along x and y for an initial source polarized along y and zo = - ,2pm. Due to a significant z 
component contribution, the y slice is slightly wider than x and the interference effects are also less pronounced. This re- 
sults in an asymmetry of the PSF in the x-y plane. There is also some asymmetry about the image in depth. This is ex- 
plained by the spherical aberration caused by the top film surface. By Snell's law , the marginal ray at 0.95 NA will cross 
the optical axis at the bottom of the film or *'* lpun. However, any paraxtaTray will cross die axis at approximately 
z' - A5[im. The difference between these represents the longitudinal spherical aberration. The fact that most of the power 
is centered around a depth of 0.5tun is probably due to the smaller amplitude of the marginal plane waves from reflection 
losses at high angles and the greater number of paraxial plane waves . 

Figure 7 shows a comparison of a the film PSF using a 0.6 and 0.95 NA system. The film is moderately absorbing with 
Tp = 1.8- .041 over a silicon substrate with t], * 1.7 - 3.381. For each NA, the film surface is offset to result in the marginal 
ray crossing the axis at a depth of Opm (no offset), 0.5pm, and 1pm. Since the substrate is no longer matched, standing 
wave interference patterns are observed. Decay of the irradiance due to absorption is also observed. The effects of 
polarization are more apparent at 0.95 NA, as evidenced by the difference in behavior between the x and y slices. 

A tri-bar object with 0.17 lines and 0.17 spaces parallel to the x axis is shown in Figure 8a. The film, substrate, and zo 
offsets are identical to the previous 0.95 NA example. The resultant y meridional slice images using a 0.95 NA system 
are presented in Figures 8b-8d for x and y polarizations. Again, interference effects are less with the y polarization due to 
the influence of the z component. In Figure 8b the onset of spurious resolution is seen at the bottom of the film by the 
presence of four major irradiance lobes instead of three. The bending nature of the iso-irradiance contours are clearly the 
result of the vectorial nature of this model. Figure 9 shows a plot of the irradiance about the zero position at a depth of 
z = 0 for the images in Figure 8b. Differences such as these become relevant in micro-lithography since the tolerances on 
the total processed structures maintained in photoresist must not exceed ± 10%. 

4. Conclusion 

A vector imaging model for thin films has been presented that is based on a plane wave decomposition of the radiation that 
propagates from the pupil, and a thin-fihn treatment of propagation into the film stack. The model is valid for NAs ap- 
proaching unity and and lends to minimize the number of plane waves needed for an image. It is shown that, when the 
illumination is polarized, asymmetries are created due to high-NA polarization effects and spherical aberration caused by 
a film surface. The image iso-irradiance contours that are presented illustrate the vectorial nature of the model 
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Fig. 1. Total absorbed power of an incident plane wave as a function of coated photoresist thickness for 
a moderately absorbing photoresist over a silicon substrate at X=436nm The angle of incidence, «, is 
given as the NA in air, i.e. sin 6 
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Imaging Lana 




Fig. 3. Coordinate system used to define the projection vector, k. 
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Fig. 4. Schematic of a plane wave propagation vector and position vectors in reference to the pupil and 
the first thin film. 
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Position (/xm) Position (/im) 

(a) <b) 

Fig. 6. Iso-irradiance contours of a point image in a \\xm photoresist film of n, = 1.8 on a matched 
substrate for a 0.95 NA system, X = 248/tm, and an initial polarization along the y axis. Meridional slices 
along (a) y axis, (b) x axis. 







Depth 0*m) Depth (/un) Depth (/im) 


l95 NA 


B \ Poeition 0*m) 

Poeition (pm) , , L l . 

Fta 7. Iso-inrad ia nce contours in of a point image in a lpm photoiesist film over a si icon su 

0.6 and 0.95 NA. The source is initially y polarized and \ = 24Snm The ' S * . ‘ ^ JJj 

along the y axis and the dotted line is a slice along the x axis. Marginal ray focus at (a) Op* depth, (b) 

0.5 pm depth, (c) 1.0pm depth. 
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Fig. 9. Pkxs of irradiance for x and y polarizations at the top surface of the film for a zero offset using 
the parameters in Figure 8b. 
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Table I 


P (a , p) 

Initial Polarization State 

Components 

X 

Y 

PxS 

P 2 

- a P 


1-7* 

1-7* 

P« 

7 a 2 

aP7 

1-7* 

1-7* 

P, 

1 * 



1 + 7 

1+7 

P,S 

-<*P 

a 1 


i-y 

1-7* 

TV 

aPv 

yP 2 

1-7* 

1-7* 

P, 

-<*P 

i-P 


1+7 

1 +7 

P* 

-a 

-P 

P* 

-a 

-p 
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